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1 Introduction 



LEPl has performed a gigantic task in testing QCD predictions. In this it has benefited from the 
very large statistics available, the substantial lack of background, and the fact that initial state 
radiation plays only a minor role on the resonance. At LEP2, QCD tests are more challenging. 
Initial state radiation is very important, there is a WW production background, and statistics 
are somewhat limited. In fig. |l]we show the annihilation cross section as a function of the centre 
of mass energy. The figure reports the Born cross section for the production of hadronic final 

10^ 
104 
10^ 

ft 
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10° 

100 200 300 400 500 

Ecm (Ecut) in GeV 

Figure 1: Hadronic cross sections as a function of the centre of mass energy. 

states through the Z/7 annihilation process, and the same cross section with the inclusion of 
the initial state radiation. This increases the cross section considerably above the Z resonance 
due to the e^e" —>■ Z 'j process, in which the hadronic system has an invariant mass equal to 
the mass of the Z boson. In the figure we also show the hadronic cross section at a fixed Ecm, as 
a function of a lower cut i^cut on the invariant mass of the hadronic system (dot-dashed line) . 
With Ecut = this cross section coincides with the value of the dotted line at 175 GeV. As 
the cut is increased above the Z mass, the cross section drops suddenly, and it approaches the 
partonic cross section at 175 GeV. As the cut approaches 175 GeV the cross section vanishes, 
but it is quite clear that if we allow for few GeV of initial state radiation, its value is very 
close to the Born cross section. Assuming therefore a 20 pb cross section, with an integrated 
luminosity of 500 pb~^ we expect 10000 hadronic events. From the figure we see that the W 
background is not a negligible one, and further cuts should be imposed to get rid of it. From 
statistics alone, the error on a measurement of the total hadronic cross section is 1%. Since 
Ctot = crtot(l + tts/TT + ...), we would expect a 25-30% error on a determination of as from the 
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hadronic cross section, not including systematics. Therefore, a useful measurement of as from 
the total cross section will not be possible at LEP2. Instead, it will be possible to determine as 
from jets. The rule of thumb in these cases is that we expect most events to be two-jet events, a 
fraction ag of three-jet events, and a fraction of 4-jet events. With 10000 hadronic events, 
we would have 1000 three-jet events, which will allow us to determine as with a statistical 
precision of 3%. Assuming that as{Mz) = 0.123, we expect as(175 GeV) = 0.112, a 10% 
variation. It seems therefore possible to see the running of as between LEPl and LEP2. 

A large fraction of this report will be dedicated to the problem of measuring as from jets at 
LEP2. In Section ^ the relevant experimental aspects of event selection and background cor- 
rections will be dealt with. In the Sections ^ and 4 the present status of theoretical calculations 
for jet shape variables will also be given. 

Using the large number of hadronic events, studies of particle spectra will certainly be 
possible. Section |^ is dedicated to fragmentation function studies at LEP2. The study of 
fragmentation functions is a relatively recent topic at LEPl. Measurements of the various 
components of the quark and gluon fragmentation functions have been performed at LEPl, 
and they allow us to make an absolute prediction for the fragmentation function at LEP2 
energies, and also for the fragmentation function in W decays. We will see that it is very 
difficult to see scaling violation effects from LEPl to LEP2. It is nevertheless important to 
measure the fragmentation function to check for the consistency of the whole approach, since 
important assumptions are often made when performing the fit (for example, flavour SU(3) 
symmetry). A study of scaling violation towards the small x region has not yet been performed 
even at LEPl, mostly because of the lack of a complete theoretical calculation. We will present 
the relevant theoretical ideas in Section |5.3| . 

Section ^ will be dedicated to the measurement of particle multiplicities at LEP2. QCD 
makes a prediction for the energy dependence of the multiplicity, and for the shape of the mul- 
tiplicity distribution, based upon the assumption known as local parton-hadron duality. The 
measurement of the multiplicity in heavy-flavoured events has recently received some attention, 
and will also be considered here. Based again upon the idea of local parton-hadron duality, 
QCD predicts many features of the small-x particle spectrum and correlations. Section |^ will 
deal with these topics. 

2 Event Selection and Event Shapes — Experimental^ 
2.1 Introduction 

A number of interesting studies of QCD may be performed at LEP2 using — > qq events. 
Although the number of events will be much smaller than at LEPl, it may be sufficient to 

^The present Section is mostly work of D. Ward, including contributions from S.Bethke, G. Cowan, D. Lanske, 
and C.Padilla. 
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explore aspects of the energy evolution of QCD. In this study we focus on the determination 
of as- The value of as(Mzo) has been determined using a number of techniques involving jet 
rates and event shape observables at LEPl|[l|, |, ||, H, ||, H and at SLD f^. For example, 
using a combination of resummed next-to-leading log (NLLA) and 0{ag) QCD calculations, 
an average measurement of 

as(Mzo) =0.123 ±0.006 



was obtained |T0|. Taking the typical centre-of-mass energy at LEP2 to be 175 GeV, we may 
expect the value of as to be reduced to 0.112. Although the change in as is not great compared 
to the uncertainty on the LEPl measurement, it should be noted that the error at LEPl is 
predominantly theoretical in origin, and thus may be largely correlated between LEPl and 
LEP2. We may therefore hope to make a useful measurement of the difference in as between 
the two energies. 

The experimental difficulties at LEP2 are somewhat different from those at LEPl. At LEPl 
hadronic Z" decays could be readily identified with efficiencies in excess of 98%, and with 
negligible background. At LEP2 there are extremely large radiative corrections, and W"^W~ 
events may contribute a significant and troublesome background. Therefore, in Sect. p.2| we 
investigate the problems of selecting a sample of non-radiative Z°/7 — > qq events, and discuss 
the extent to which these selection procedures may bias the events selected. 

It will turn out that the events which may be selected most cleanly are those nearer to 
the two-jet region. Multi-jet events are much more susceptible to contamination from W"'"W~ 
events. Since statistics are also meagre, and most of the events lie in the two jet region, this 
suggests that techniques based on the resummed NLLA QCD calculations will be most effective 
in determining as, since these calculations are expected to describe the two-jet region best. We 
have therefore focused on those event shape variables for which complete resummed NLLA 
calculations are available, namely Thrust (T), heavy jet mass [Mh), total jet broadening (Bt) 
and wide jet broadening {B\y) [11]. We also examine jet rates in the Durham jet-finding scheme, 
for which NLLA calculations are available - specifically the observable 1/2^^ which is the value of 
ycut at which the event changes from two- to three-jet. All these variables are discussed in, for 



example, Ref. 0. The next-to-leading order calculation for ^23''' as used so far by the LEP 
experiments was known to be incomplete. Recently, however, a more complete calculation has 
been presented [0. We have not yet studied this new calculation, for compatibility with the 
existing LEPl results. Resummed calculations are also now available for the C-parameter [P 
though they are not yet published, and are therefore not discussed here. 



2.2 Selection of — ^ qq events 



The discussion here will be based on events generated with Pythia [|T^] version 5.715, with 



hadronization parameters tuned to LEPl data [jT6l. The examples given below will relate to 
events processed through the Opal detector simulation, but it is to be expected that similar 
results would hold for the other experiments. The cross-sections predicted for Z°/7 qq events 
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and for the principal source of background W+W — ^-qqqq are as follows: 



Reaction 


Cro 
161 GeV 


ss-section / 
175 GeV 


pb 

192 GeV 


e+e~ W+W" qqqq 


1.69 


6.34 


7.92 


e+e~ Z°/7 — > qq 

e+e- Z77 qq ; Ei^r < 30 GeV 

e+e- Z77 qq ; E,,r < 1 GeV 


149.6 
39.0 
26.1 


116.8 
29.9 
20.0 


90.6 
22.4 
15.1 



The Z°/7 — s> qq cross-section is also given for two cuts on the amount of energy lost in initial 
state radiation. The useful cross-section for QCD studies is the non-radiative cross-section. The 
cut at 30 GeV corresponds roughly to the minimum in the hadronic mass spectrum da/dMh 
between the non-radiative process and radiation down to the Z° pole. Unless otherwise stated, 
the results shown relate to 175 GeV. 

It is helpful to consider the selection of Z°/7 — * qq events in two stages. In stage I we 
remove the leptonic and highly radiative events, and the W+W-— i>qq£r^ events, mainly using 
cuts on multiplicity and energy /momentum balance. These cuts introduce rather little bias 
into the Z°/7 — > qq event sample. The stage II cuts are to remove W+W-^qqqq events, and 
are more problematic, since they turn out to bias the selected Z°/ 7 — > qq sample significantly. 

Typical stage I cuts would be as follows: 

• Require |cos^t| < 0.9 to ensure reasonable containment of the event, where 6t is the 
polar angle of the thrust axis. 

• Require the number of charged tracks to be N^h > 6 to remove purely leptonic events. 
This cut causes a negligible loss of Z°/ 7 ^ qq events. 

• In Fig. 1^ we plot Rms against Rmiss for various classes of events, where Ryis is the visible 
energy scaled by the centre of mass energy Ecm., and Rmiss is the missing momentum 
scaled by Ecm.- It is desirable to have the best resolution on the visible energy and 
missing momentum, which involves using an algorithm to combine the information from 
the charged tracks, electromagnetic and hadronic calorimeters so as to reduce double 
counting. We note that the non-radiative Z''/7 —>■ qq events and the W+W-— >qqqq events 
are peaked around Rms = 1 and Rmiss = 0. The W+W-— >qq£r^ events, and most of the 
radiative Z''/7 qq events lie away from this point. Typical cuts are shown by the lines 
in Fig. |[ 

• Fig. ^(b) reveals a group of radiative Z°/7 — > qq events having R^s ~ 1 and Rmiss ~ 0. 
In these events, the radiative photons are detected in the electromagnetic calorimeter. 
Such photons may be identified using standard criteria on lateral shower shapes. The 
cluster should also be required to be be isolated, for example by demanding that within 
a cone of half angle 0.2 rad centred about the cluster less than 1 GeV is observed. If 
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Figure 2: Plots of Rms against Rmiss for (a) Z^/j — > qq events having less than 30 GeV initial 
state radiation (b) Z°/7 qq events having more than 30 GeV initial state radiation (c) 
W+W~— >qqqq events (d) W+W~— >qq£l^ events. The hnes show typical cuts. These plots are 
at Ecm — 175 GeV, though they are only weakly energy dependent. 
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the energy of the most energetic cluster satisfying the above criteria exceeds 0.6 x p^, the 
event is rejected. Here, is the expected photon momentum in an e"'"e~ —>■ Z°7 event, 
i.e. p, = iEl^,-Ml)/2E,,^,. 

The cross-sections for the various channels of interest before and after these stage I selection 
cuts are listed in table |I|. Hence, at Ecm = 175 GeV, in the region | cos6't| < 0.9, the stage I 



Channel 


Cross-section /pb 
cos 6't < 0.9 


Cross-section /pb 
after stage I cuts 


Z77 qq (Eisr < 1 GeV) 


17.90 


16.51 


ZVt qq {Eisr < 30 GeV) 


26.68 


23.99 


Z77 ^ qq {Eisr > 30 GeV) 


73.37 


1.07 


W+W-^qq£z7^ 


5.88 


0.06 


W+W"^qqqq 


6.08 


5.74 



Table 1: Cross-sections at Ecm = 175 GeV, based on Pythia. 



cuts accept 92% of the non-radiative Z°/7 — > qq events, whilst accepting only around 1.5% 
of the radiative qq events and W^W^ ^qqiV^ events. The W"*'W~^qqqq events are 

accepted with high efficiency. The corresponding figures at 192 GeV and 161 GeV are essen- 
tially the same. Backgrounds from two-photon events, Ze"'"e~ and Weu final states have been 
examined, and appear to be negligible. ZZ— i>qqqq does contribute, but at a much lower rate 
than W"'"W'^qqqq, with similar characteristics. 

The main feature which distinguishes the W"'"W~— s>qqqq events (and also the much smaller 
contribution from ZZ— >qqqq) from the Z°/7 —>■ qq events is that the former contain four quarks, 
and thus generally have four or more jets, and are therefore less coUimated. Furthermore, the 
invariant masses of appropriate pairs of jets should equal the mass of the W boson. We have 
examined the use of the following variables in separating these event classes: 

• The "narrow jet broadening", B]\f. The event is divided into two hemispheres, S±, by 
the plane orthogonal to the thrust axis, tit- In each hemisphere, the quantity B± = 
Z]ie5± \Pi ^ '^T|/2I]j \Pi\ is computed, where the sum in the denominator runs over all 
particles, whilst that in the numerator runs over one hemisphere. i?7v is defined by 
Bn = min(5+,5_). 

• The scaled "light hemisphere mass". Ml/ Ems- The event is divided into two hemispheres, 
S±, by the plane orthogonal to the thrust axis, and the invariant mass of each is computed, 
M±. Then, Ml is defined by Ml = min(M+, M_). 

• The value of ycut at which the event changes from 3-jet to 4-jet in the Durham jet finding 
scheme, yi^^ . 
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• Using the Durham jet finder, the event may be forcibly reconstructed as having four 
jets. The invariant masses of pairs of jets may be formed, from which we define the 
variable D = min [{Mij -Mwf + {Mki -Mwf] where the minimum is taken over the 
permutations {ij; kl) = (12; 34), (13; 24), (14; 23). Various ways of scaling the jet energies 
in order to improve the W mass resolution have been considered in connection with the 
W mass determination, but have not been used here. 

In Fig. 11(a) we show the distributions of for non-radiative Z°/7 — > qq events {Eisr < 
1 GeV) and for W"'"W~^qqqq events, after the stage I cuts. In order to judge the correlation 




l-T -I 

0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 



Figure 3: (a) Distributions (at 175 GeV) of Bj^ for Z^/j qq events having less than 1 GeV 
initial state radiation (open histogram) and for W+W"— >qqqq events (shaded) (b) average 
values of (1 — T), Bw and 1/23 ^ (scaled by their overall mean values) as a function of Bn 

between B^ and the observables which we would wish to use for the determination of Og we 
show in Fig. ^(b) the average values of (1 — T), B^ and 1/2^'' (normalized to their overall 
mean values) for non-radiative Z^/j — >■ qq events as a function of B^. It is evident that the 
W+W^^qqqq contribution can be reduced to almost any level desired by cutting on B^, but 
at an increasing cost in bias, and a corresponding loss in statistics. Generally, Ml and 
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show similar behaviour to B^- The variable offers a less clean separation between the 
Z°/7 qq and W~'"W~— >qqqq events, but it appears that it may introduce somewhat less, or 
different, bias, and may thus be complementary. 

These observations may be quantified in table ^ below, where we show the effect of various 
possible stage II cuts on the Z°/7 qq non-radiative signal and the W"''W~^qqqq background. 
We give the average values of 1 — T, Bw and as an indication of the bias caused by the 
cuts. We note that the stage I cuts cause only a small bias. We show several possible cuts 





Z77 ^ qq{E,sr < 1 GeV) 


W+W ^qqqq 


Cut(s) 


(T /pb 


< i-r> 


< Bw > 


< y^.^^ > 


0" /pb 


cos 6^2^ < 0.9 


17.90 


0.0598 


0.0708 


0.0196 


6.08 


Stage I only 


16.51 


0.0587 


0.0701 


0.0195 


5.74 


Bn <0.07 


15.73 


0.0528 


0.0668 


0.0169 


1.34 


Bn <0.06 


15.27 


0.0504 


0.0653 


0.0160 


0.87 


Bn <0.05 


14.53 


0.0474 


0.0632 


0.0149 


0.50 


Bn <0.04 


13.26 


0.0433 


0.0601 


0.0133 


0.24 


Ml <0.175 


15.26 


0.0504 


0.0656 


0.0161 


1.37 


<0.0065 


15.34 


0.0501 


0.0644 


0.0156 


0.88 


^2 >300 GeV2 


14.17 


0.0543 


0.0667 


0.0164 


1.90 


^2 >600 GeV2 


12.25 


0.0504 


0.0638 


0.0140 


0.89 


L>2 >600 GeV^ and Bn <0.06 


11.71 


0.0460 


0.0613 


0.0124 


0.24 


D'^ >300 GeV2 and Bn <0.05 


12.79 


0.0457 


0.0615 


0.0132 


0.24 


Bn - 7/^/2000 < 0.03 


13.84 


0.0451 


0.0612 


0.0137 


0.25 



Table 2: Cross-sections at 175 GeV accepted after the Stage I cuts, and after various possible 
Stage II cuts. The average values of various relevant observables are also shown, to indicate 
the level of bias introduced. 



on the Bn variable. The background from W^W^— >qqqq may be reduced, for example, to a 
level of 4% with an efficiency for selection Z°/7 — > qq events of 82%. However, the sample of 
Z°/7 qq events accepted is strongly biased. The bias, as measured by the mean value of the 
observable, tends to be greatest for y2^^ and smallest for B^y. We show similar results for cuts 
on Ml and 2/3^^ where we have chosen cuts which yield roughly the same Z°/7 qq efficiency 
as the Bn < 0.06 cut. Cutting on Ml is less effective than Bn at removing W"'"W~— >qqqq 
background, while a cut on y^^^ gives essentially the same performance as -Bat. The cut on 
> 600 GeV^ yields the same W"'"W~^qqqq contamination (7%) as the Bn < 0.06 cut, 
but for a significantly lower Z°/7 qq efficiency (69% compared to 85%). Using yields a 
somewhat smaller bias on 1 — T, but the bias on y2^^ is a little greater. The two observables Bn 
and are not strongly correlated (whereas, for example, Bn and Ml are highly correlated), 
suggesting that a joint cut on the two variables could give better separation. Examples are 
given in Table |^. The W+W"— >qqqq background may, for example, be reduced to around the 
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2% level for a qq efficiency of almost 80%, with somewhat less bias than a cut on Bn 

alone. The precise cuts chosen for the separation of Z°/7 —>■ qq and W~^W^^qqqq events may 
therefore need to depend on the analysis being performed - whether a high purity is demanded, 
or whether a comparatively unbiased sample is required. 

In Fig. |^(a) we show the distributions of a typical observable which may be used for the 
determination of as, (1 — T), after the stage I cuts. We compare the Z°/7 — > qq non-radiative 
{Eisr < 1 GeV) signal with the W+W~— >qqqq background. In Fig. ||(b) we show the same 
distributions after the stage II cuts, taking < 0.05 as a typical stage II cut. As expected. 




Figure 4: (a) Distributions (at 175 GeV) of (1 — T) after the stage I cuts. Z°/7 qq non- 
radiative {Eisr < 1 GeV) events are shown as points with errors, and W+W"— s>qqqq events by 
the shaded histogram, (b) as (a), after applying the stage I cuts and the stage II cut B]\r < 0.05. 
(c) Biases to the distribution of (1 — T). The stage I cuts and the stage II cut B^ < 0.05 are 
applied. The closed points show the fraction of Z°/ 7 — qq non-radiative {Eisr < 1 GeV) events 
accepted after cuts. The open points show the ratio of all accepted Z°/7 qq events after cuts 
to non-radiative {Eisr < 1 GeV) Z°/7 — > qq events before cuts. 

the background tends to be concentrated toward large values of (1 — T), i.e. the region of 
hard gluon emission in the ZP/j qq reaction. The two-jet region of the Z°/7 — >• qq process 
is relatively free of background. Other stage II cuts give similar results. In Fig. H(c) we show 
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the efficiency of the stage I+II cuts, taking Bj^ < 0.05 as a typical stage II cut, as a function 
of the (1 — T), for Z°/7 qq non-radiative {Eisr < 1 GeV) events (sohd points). As expected, 
the cuts bias against large values of (1 — T). In Fig. H(c) we also show as open points the ratio 
of the distributions of all accepted Z^/'y — > qq events (including radiative events) to those of 
the non-radiative events before selection cuts. In general the effect of initial state radiation is 
to bias the distribution towards higher values, but this is counteracted by the tendency of the 
cuts to reject events with high values of the observables. The net effect is that the distribution 
of the accepted radiative Z°/7 — > qq events is quite similar to the distribution of non-radiative 
events before cuts, and so the ratios in Fig. ^(c) are increased roughly uniformly. Other stage 
II cuts give similar results, though the efficiencies may be systematically higher or lower. 



2.3 Determination of cks 

Before comparing with QCD calculations, the observed data must be corrected for the effects 
of detector resolution, the acceptance of selection cuts and the effects of background (Fig. 
The influence of hadronization must then be accounted for, and one standard way of doing this 
is to multiply the corrected hadron level data by the ratio of the parton level to hadron level 
distributions from a Monte Carlo model. In Fig. ^ we show these ratios for (1 — T), based on 
Jetset7.4, at 175 GeV (LEP2) and 91.2 GeV (LEPl). We note that the correction factors 




Figure 5: Hadronization corrections for the distributions of (1 — T). 

at LEP2 are significantly closer to unity, especially at small values of (1 — T), corresponding 
to the two-jet region. Similar comments apply to the other observables. A requirement for a 
credible analysis is that the correction factors be not too far from unity. 

In this study, we investigate three types of QCD calculations, which may be used as the 
basis of a measurement of as from event shape variables. These are: 



0(ag ) The QCD matrix elements, expanded as a power series in ag are fully known to O^a"^) |lj 
From previous studies at LEPl we know that these calculations are applicable in the "3- 
jet" region, i.e. the region dominated by hard gluon radiation. A significant uncertainty 
in applying the Ola"^) calculations is the choice of renormalization scale, /i, represented by 
Xfj, = fi/Ecm.- The region over which the data can successfully be ffited can be extended 
further into the 2-jet region by choosing a small value of ~ 0.1 (typically). 
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NLLA In the 2-jet region, the expansion in powers of as is bound to fail, because large 
logarithms arise associated with coUinear and soft gluon emission. In this region, "NLLA" 
calculations are available which resum the leading and next to leading logarithms to all 
orders in as- It has been shown in ref. that such calculations may be used to derive as 
at LEPl, but that it is necessary also to include a sub-leading term of the form G2ias^L 
in order to achieve a good description of the data. 

Combined NLLA-|-(9(ag) The most complete embodiment of our present knowledge of QCD 
comes from combining the 0{a1) and NLLA calculations. It is necessary to match the 
calculations in such a way as to eliminate double counting of terms, and there are several 
ways of doing this. These have been studied at LEPl, based on which we choose the 
"Ini?" matching scheme for the present work. 

To assess the range of validity of these calculations at LEP2 we proceed in the following 
empirical manner. We have generated distributions of the five observables, (1 — T), M^^, Bt-, 
Bw and y^^\ at the parton level, using the Jetset7.4 parton shower model without initial state 
radiation. We can assume that the data, after correction for detector acceptance, the effect of 
selection cuts and background, and hadronization, would closely resemble these distributions. 
For each observable, we then determine the largest range for which the theoretical calculations 
reproduce those from jETSETwith an acceptable x^/DOF. The results are summarised in 
Table ^. We note that the Oi^al) calculations may (in most cases) be extended to lower values 
of the observables by fitting x^. The NLLA or combined calculations allow a description down 
to still lower values, but, particularly in the case of the pure NLLA calculations, the higher 
values of the observables are less well modelled. The NLLA and combined calculations for Bw 
tend to give a rather poor description of the Jetset "data" (as seen at LEPl). The pure NLLA 
calculations are not applied to y'^\ since they are known to be incomplete, and in fact yield a 
poor fit to the Jetset distributions. 



Observable 


Oia!) ix,=l) 


0{al) {x^ fitted) 


pure NLLA 


Combined 0(a2)+NLLA 


(1-T) 


0.09-0.3 


0.05-0.3 


0.02-0.17 


0.02-0.3 


Mh 


0.20-0.55 


0.14-0.55 


0.10-0.35 


0.14-0.55 


Bq^ 


0.11-0.3 


0.10-0.3 


0.05-0.18 


0.05-0.22 


Bw 


0.06-0.26 


0.06-0.26 


0.02-0.12 


0.05-0.17 


(D) 
1/23 


0.015-0.2 


0.005-0.2 




0.005-0.2 



Table 3: Approximate ranges of applicability of various types of QCD calculation. 



If, for example, we require that the hadronization corrections lie between 0.8 and 1.2, that 
the Z°/7 — i> qq acceptance be greater than 50% and that the W"'"W~— s>qqqq contamination 
be less than 50%, the regions where the data can be used reliably would be roughly 0.03-0.2 
for (1 - T), 0.15-0.4 for Mh, 0.06-0.2 for Bt, 0.03-0.18 for Bw and 0.005-0.09 for yl^^K By 
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comparison with Table ^ it is evident that the regions in which rehable data may be obtained 
are best matched by the regions in which the combined NLLA+(9(ag) calculations are valid. 
Since these are also the most complete calculations, this would appear to be the most promising 
approach. 

We next assess the precision on ag which could be achieved using 500 pb~^ of data at 
LEP2. In order to do this, we take the Jetset7.4 parton level distribution, with statistical 
errors corresponding to this integrated luminosity (approximately 6500 Z°/7 —>■ qq events). 
We then fit the QCD theory to infer as, fitting in the range of the observable given by the 
overlap of the ranges in Tables ^ and the regions where reliable data may be obtained. A 
typical fit (of the 0{a'^)+NLLA calculations to (1 — T)) is shown in Fig. |^. We find that 
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Figure 6: Typical fit of the(9(a;s)+NLLA QCD calculations to (1 — T) in order to determine 
as- The dotted lines delimit the fit region. 

the 0{a'^) calculations yield typical statistical errors of ±0.0024, which are larger than the 
NLLA and combined 0{a'^)+NLLA calculations (typically ±0.0015) because the former are 
only applicable towards the 3-jet region, where the few events are found. It also appears that 
the statistics are generally insufficient to permit a precise determination of the scale factor 
for the 0{a'^) fits. The pure NLLA and combined 0(ag)±NLLA calculations both appear to 
be competitive, and offer the possibility of measuring ctg with a statistical precision of around 
±0.0015. For the event shapes (1 — T), Mh, Bt and Bw, the NLLA tend to yield smaller 
values of as, and the 0{a'^) calculations larger values; the same trend was noted at LEPl |Q. 

As at LEPl, the combined (9(ag)±NLLA method will probably be the preferred technique, 
because it represents the most complete theoretical calculations, and allows the largest fraction 



14 



of the data to be included in the analysis. For the discussion of possible systematic uncertainties, 
we therefore focus on these calculations. In ref. 0, for example, a wide range of systematic 
effects were investigated. The largest contribution was found to arise from variation of the 
renormalization scale factor x^. Other significant effects arose from varying the hadronization 
model, particularly from the use of the Herwig model, and from the influence of b-quark mass 
effects. We have estimated the systematic errors resulting from these effects at 175 GeV, and 
compared with the LEPl experimental results. 



• The renormalization scale factor is varied in the range 0.5< <2.0. The changes in as 
are highly correlated with those at LEPl, though about 20% smaller on average. If we 
assume that it makes sense to choose the same scale factor at LEP2 as at LEPl, then 
the effective systematic uncertainty on the change in as between LEPl and LEP2 would 
be about ±0.0015 for (1 - T), ±0.0025 for Mh, ±0.0015 for Bt, ±0.0005 for Bw and 
±0.0003 for ySl 

• The influence of b-quark mass effects may be crudely accounted for by basing the parton 
level distributions in the correction procedure only on udsc quark events. At LEPl this 
correction was found to increase as by about 0.002 for most observables. Not surprisingly, 
the effect is much smaller at LEP2. However, the relevant point is the difference between 
the LEPl and LEP2 uncertainties, which is of the order of 0.002 (somewhat larger for Bj- 
and smaller for M^)- 

• The Herwig model offers a quite different hadronization scheme from Jetset. Since 
the hadronization corrections are smaller at LEP2 than at LEPl, we would expect the 
uncertainty associated with the use of different models to be reduced. This is generally 
the case, but the correlation between the Herwig uncertainties at LEPl and LEP2 is 
unclear. This is partly because different fit regions have been used, and also different 
versions of the models. Clearly, in order to establish a reliable systematic uncertainty on 
the difference in as between LEPl and LEP2 it would be necessary to make a more careful 
analysis using consistent versions of the models at the two energies. For some observables 
at least (e.g. Bw and yis ^) it seems plausible that the hadronization uncertainty could 
be quite small. 



In summary, it appears that systematic errors would not preclude making a useful measure- 
ment of the difference in as between LEPl and LEP2. The renormalization scale uncertainty 
seems to be comparable with or smaller than the statistical error. The uncertainty associated 
with b-quark mass effects could perhaps be reduced by further analysis and theoretical work. 
The uncertainties associated with the choice of hadronization models are less clear; it may 
be necessary to reanalyse the LEPl data using the same models and parameter sets as em- 
ployed in the LEP2 analysis, and the same fit regions, in order to minimise the uncertainties. 
Nonetheless, it seems that the systematic errors could be quite small, for some observables at 
least (especially Bw and y^sK according to our study). It may be noted that recent studies of 
non-perturbative (power) corrections to the mean values of event shape observables suggest 
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that certain observables or combinations of observables might be expected theoretically to have 



especially small hadronization uncertainties (e.g. j/23 or T — 2C/3n [T^ 



3 Event Shapes — Theoretical^ 



Since the completion of the Yellow Report for LEPl much progress has been achieved 
in the theoretical calculations of shape variable distributions. A technique of resummation of 
contributions enhanced near the two-jet region has been studied and fully implemented in refs. 



m 12, 



mm 



Furthermore, new calculations of shape variable 



distributions (implemented as computer code) have become available. 

Calculation of shape variables are all based upon the original work of ref. [17]. This calcu- 
lation was also performed in ref. |2T|. Although the analytic results did agree, several problems 
where found in the comparison of numerical results (see ref. |2y] for a small review). While 
at the time of ref. it was hard to find precision calculations of jet shape distributions 
that agreed with each other, today we have at least three general purpose programs that do 
agree. One, the program EVENT, was developed for ref. 
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Results of shape variables dis- 
tributions performed with this program are reported there, and have served as a benchmark 
for comparison with other computations. In ref. |^ a new computation was performed, which 
agrees with good accuracy with ref. Furthermore, very recently, yet another calculation 
was completed In ref. also oriented events are implemented, and apparently they will 



also be implemented in ref. This means that it will be possible to compute distribution 
of shape variables that do depend upon the orientation of the incoming beams axis, unlike all 
shape variables that where used up to now (see the next section). 

The most disturbing disagreement on shape variables was found to be on the Energy-Energy 
correlation (EEC). The computation performed in ref. |Q was found in important disagreement 
with other calculations, and in particular with ref. |20]. Recently, in ref. [35| the calculation 



of ref. pi] was repeated. The result of the new calculation was found in disagreement both 

|20[| and with ref. No clear statement is made in ref. upon the 



with the result of ref. 



1 with ref. No clear statement is made in ref. 
origin of the discrepancy. It is however claimed that the disagreement comes from the region 
in which besides the quark-antiquark pair, two soft gluons have been radiated. The EEC is 
in fact peculiar, in the sense that even configurations with thrust near 1 can contribute to the 
EEC at angles far away from and vr. Because of the lack of a more complete theoretical 
paper form the authors of ref. ||35|, we thought that the most useful thing to be done for the 



present report is to perform a high-precision comparison of the different computations of the 
EEC, that can serve as benchmark for future calculations. In order to achieve high precision, 
instead of computing the EEC itself as a function of the angle, we computed its moments. The 



^Written by P. Nason, including contributions of M.H. Seymour, N. Glover and K. Clay. 
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energy-energy correlation is defined as 



We define 



/ EEC(x) sin^^- X cos- xdcosx= ^4e'? + (g) ' ^Eb"^ + Oia^') (2) 

where = Q;s(-Ecm)- The coefficients -Beec^ have the following colour structure 

4ec"^ = Cf {Ca 4:'"^ + Cf + Tf nj 47'")) . (3) 

We then asked K. Clay (C), N. Glover (G), M. Seymour (S), and the author (N), to compute 



^(m,n)^ ^(m,n)^ j^irn,n) ^j^]^ ^^iq programs of rcf. , , and for m = 0, . . . , 5 and 



n = 0,1. All four computations agreed within errors for the Bj^'^^ term. In the other two cases 
we found disagreements. The results are reported in tables ^ and ^ 



m 


n 


N 


G 


S 


C 








50.82 ±0.05 


50.54 ±0.03 


50.72 ±0.02 


46.4 ±0.2 


1 





35.76 ±0.04 


35.53 ±0.02 


35.64 ±0.02 


32.09 ±0.06 


2 





28.94 ±0.03 


28.75 ±0.02 


28.82 ±0.02 


25.73 ± 0.04 


3 





24.92 ±0.03 


24.75 ± 0.02 


24.80 ±0.02 


22.03 ± 0.04 


4 





22.20 ±0.03 


22.05 ±0.02 


22.09 ±0.02 


19.54 ±0.04 


5 





20.21 ±0.03 


20.07 ±0.02 


20.10 ±0.02 


17.74 ±0.03 





1 


-6.468 ±0.006 


-6.50 ±0.01 


-6.455 ± 0.005 


-6.0 ±0.15 


1 


1 


-2.356 ±0.004 


-2.365 ±0.009 


-2.344 ±0.003 


-2.15 ±0.03 


2 


1 


-1.189 ±0.003 


-1.194 ±0.008 


-1.177 ±0.003 


-1.06 ±0.02 


3 


1 


-0.714 ±0.003 


-0.718 ±0.007 


-0.702 ±0.003 


-0.62 ±0.01 


4 


1 


-0.478 ± 0.003 


-0.479 ±0.007 


-0.466 ± 0.003 


-0.41 ±0.01 


5 


1 


-0.344 ±0.003 


-0.344 ±0.006 


-0.331 ±0.003 


-0.28 ±0.01 



Table 4: Comparison of different computations of the B^ ' coefficients. 



It is clear that the results N, G and S agree with each other with high accuracy, while C 
is seriously different. Observe that, although for all practical purposes N, G and S agree with 
each other, there are among them discrepancies of several standard deviations. We attributed 
these differences as an underestimate of the errors, rather than to a real difference in the 
calculation. More details on the different characteristics of the three computer codes are given 



in the generator's section [38 
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-1.65 ±0.01 


-1.81 ±0.03 


-1.67 ±0.03 


8.22 ±0.02 


5 





-1.36 ±0.01 


-1.51 ±0.03 


-1.39 ±0.03 


7.69 ±0.02 





1 


4.906 ±0.002 


4.92 ±0.01 


4.892 ±0.006 


2.6 ±0.2 


1 


1 


0.240 ±0.002 


0.259 ±0.008 


0.232 ±0.004 


-0.58 ±0.02 


2 


1 


-0.383 ±0.002 


-0.367 ±0.006 


-0.386 ±0.004 


-0.80 ±0.01 


3 


1 


-0.458 ±0.002 


-0.445 ± 0.005 


-0.459 ± 0.003 


-0.72 ±0.01 


4 


1 


-0.428 ±0.002 


-0.417 ±0.005 


-0.429 ± 0.003 


-0.606 ±0.007 


5 


1 


-0.381 ±0.002 


-0.371 ±0.004 


-0.381 ±0.003 


-0.510 ±0.006 



Table 5: Comparison of different computations of the B^'^^ coefficients. 

4 Next-to-leading Order Calculations of Oriented Event 
Shapes^ 

At LEP2, it will become increasingly important to be able to cut out some angular regions to 
control the backgrounds, and to define event shapes that are invariant under boosts along the 
beam direction to study continuum events with initial-state radiation. To make predictions for 
such quantities it is essential to use the full matrix elements for e+e" qqg, including the full 
Z/7 interference and the polarisation of the exchanged boson. Two programs have recently 
become available that include these matrix elements, EERAD||39|| and EVENT2[37]. These use 
completely different methods to implement the cancellation of poles between real and virtual 
contributions, as described in 153], but the results are in excellent agreement with each other. 



As an example of an oriented event shape we study the thrust distribution as a function of 
the thrust axis direction. As usual [20], we parametrize the distribution as 



1(1 - T) ^ f , = ^A{T, cose) + 



2 r 



,2 



A{T, cos ^)27r6o log ^ ± B{T, cos I 
s 



(To dTdcosO 2tc \ 27r _ 

(4) 

The definition of thrust has a forward-backward ambiguity, so we are at liberty to define 
cos 6' > 0. The leading order term is known analytically||40||. 



A(T,cos^) = 



2(3r^-3r±2) , 2T-1 T^^n^ 2 
f 1-T ~ ^ - 2)(2 - T) ^ 1(1 ± cos I 



Author: M.H. Seymour 
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Figure 7: The coefficients of the thrust distribution for five bins in cos 6', where 9 is the angle 
between the thrust axis and the beam. The errors shown are purely statistical and are similar 
for each histogram, so we only show them for one. 



+ 



2(3r-2)(2-T)(l -T) 
2^2 



1 - 3 cos^ 9) 



(5) 



Integration over cos^ immediately gives us the expression in [20]. Numerical results for A and 
B are shown in Fig. ^ Combining these coefficients with an value and factorisation scale 
choice, = s) = 0.120, we obtain the predictions shown in Fig. ^a. Alternatively, we can 

divide out the trivial dependence on cos 6' by normalising each curve to the number of events 
in that bin, given by|50] 



1 da 
(To rfcos 9 



f(l + cos2| 



1 + 



+ 



vr 



Slog- 



1 — 3 cos 



cos 



a.. 



TT 



(6) 



where the approximation is good to better than 1%. The result is shown in Fig. §0, where we 
see that the majority of the cos 6' dependence in Fig. ^ was from this dependence of the total 
event rate and the residual dependence is rather small. Nevertheless, it should be measurable 
with the full statistics of LEPl. 



5 Fragmentation functions^ 

The measurement of fragmentation functions at different energies and the comparison with 
the theoretical predictions, either implemented in the Monte Carlo programs or deduced from 
other measured data, can be used to perform different QCD tests and to tune the parameters 
describing the fragmentation processes inside the Monte Carlo programs. 

^5 Author: C. Padilla. 
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Figure 8: Predictions for the thrust distribution for five bins in cos^ normahzed to (a) the total 
number of events and (b) the number of events in each bin. 

At the energies available at LEP II the scaled energy {x = 2Ej^^ distributions for charged 
particles can be measured for qq events in which the mass of the hadronic system is close to 
the centre-of-mass energy of the collision. Furthermore, the fragmentation function of the W 
boson can also be measured and compared to the expectation that comes from the measurement 
of the fragmentation functions for different enriched flavour samples at LEP I, after correcting 
for the small scaling produced for the different masses of the Z and the W boson and for the 
different flavour composition. 

This Section describes how the measurement of the scaled energy distributions can be made 
and what can be expected in the measurement of from scaling violations. 



5.1 Measurement of scaled energy distributions 

The measurement of the charged scaled energy distributions will follow the same procedure 
used at lower centre-of-mass energies. At centre-of-mass energies of the Z mass, hadronic 
events can be selected with very high purity and small backgrounds (coming mainly from r 
events). At LEP II centre-of-mass energies, most of the qq events (more than 75%) will radiate 
an initial state hard photon such a way that the effective centre-of-mass energy of the collision 
will be reduced to below 120 GeV. These events have a high boost along the collision axis and 
have to be removed. 

The selection of hadronic events will follow a procedure very similar to the one presented in 



subsection p.2i After some minimal requirements on track quality, number of tracks and total 
measured energy of these tracks, additional selection variables have to be considered. Good 
containment of the events can be obtained with cuts in the sphericity or thrust axis. 

Monte Carlo simulations performed in ALEPH, based upon DYMU3 and JETSET, including 
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full simulation of the detector response, show that requiring a visible mass of the event above 
120 GeV and a normalized balanced momentum of the charged tracks along the beam axis 
below 0.3, a selection efficiency of ~ 18% can be achieved. The percentage of selected events 
such that the invariant mass of the propagator is below 120 GeV is reduced to approximately 
7% with this selection procedure. 

The backgrounds from dilepton events are small at this level. However, the background 
from WW events could still be substantial. A cut in missing momentum will remove most of 
the events in which one of the W has decayed leptonically. The remaining events in which 
both W decay hadronically can be removed by considering appropriate shape variables. Events 
resulting from the fragmentation of two W bosons will have a four-jet topology that makes 
them more spherical than the ones resulting from Z/7 — qq. In subsection \i.2\ a discussion 
of the various possible approaches is given. For the case of the measurement of the scaled 
energy distributions, a cut on thrust T > 0.925 would be appropriate, since (unlike the case of 
shape variables) such a cut does not introduce strong biases in the shape of the fragmentation 
function. 

The whole selection procedure should result in a cross section for qq events of ~ 11 pb 
with less than 1% of events with the effective centre-of-mass energy below 120 GeV and with 
a background of WW events below 5%. Assuming an integrated luminosity of 500 pb~^ the 
expectation is to have ~ 6000 selected hadronic events. 

It can be assumed that the background can be subtracted statistically using Monte Carlo 
techniques, and that the distribution is corrected using a hadronic event generator (with param- 
eters adjusted to describe the data) for the effects of geometrical acceptance, detector efficiency 
and resolution, decays of long-lived particles (with r > 1 ns), secondary interactions and resid- 
ual initial state photon radiation. The bin-to-bin correction factors are below 10% using the 
selection described above. 

Figure ^ shows the Monte Carlo scaled energy distribution for the statistics of 6000 events. 
The energies of the particles before detector effects have been used to construct the distribution. 
Additional systematic uncertainties coming from possible discrepancies between the real detec- 
tor performance and the simulated one and from the dependence on the hadron production 
model used to correct the data for detector effects will have to be considered. 

The measurement of the W fragmentation function will require the selection of hadronic W 
events. The events in which one of the W decays leptonically can be selected using missing 
momentum or tagging a high-momentum lepton. The rest of the particles can be used to 
determine the momentum of the hadronically decaying W boson and to construct the scaled 
energy distribution, after boosting the particles into the rest frame of the parent W boson. In 
the case that both W bosons decay hadronically the techniques used in the measurement of the 
W mass can be used to unambiguously assign the jets to the corresponding W bosons. 
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Figure 9: Scaled energy distribution for hadronic events at a centre-of-mass energy of 180 GeV 
according to the JETSET Monte Carlo. The error bars correspond to the statistics of 6000 
events. 

5.2 Scaling violations: QCD tests 



The analysis of scaling violations with the data available at LEP and data from lower centre- 
of-mass energy experiments (PEP, PETRA, TRISTAN) has focused on the measurement of 
The prediction of scaling violations in fragmentation functions of quarks and gluons 



41, 42 



a 

is similar to that predicted in structure functions in deep-inelastic lepton-nucleon scattering. 

In an electron-positron collider, scaling violations are observed in the dependence of the 
distribution of the scaled energy of final-state particles in hadronic events on the centre-of- 
mass energy ^/s. This comes about because with increasing y/s more phase space for gluon 
radiation and thus for final-state particle production becomes available, leading to a softer 
x-distribution. As the probability for gluon radiation is proportional to the strong coupling 
constant, a measurement of the scaled-energy distributions at different centre-of-mass energies 
compared to the QCD prediction allows one to determine the only free parameter of QCD, as- 
A recent review of the relevant theoretical ideas has been given in ref . . For another recent 



theoretical analysis see [44 



A reliable measurement of scaling violations has to disentangle the true QCD evolution 
from effects due to the dependence of the flavour composition upon the centre-of-mass energy. 
Since heavy flavours, after their decay into light particles, typically have softer fragmentation 
functions, when going from centre-of-mass energies below the Z mass towards the Z mass, the 
h content increases, and it decreases again when going towards higher energies. To analyse 
the data in a model independent way, final-state flavour identification and a measurement of 
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the gluon fragmentation function are needed. This procedure has been followed in the analysis 



performed in ref. where enriched uds-, c-, and 6-quark scaled energy distributions, together 
with the measurement of the gluon fragmentation functions and the longitudinal cross section, 
have been used to constraint the fragmentation functions for the different flavours and the 
gluon. It was assumed that the fragmentation functions of u, d, and s quarks are the same. In 
the analysis presented there, a total of 15 parameters besides a.^ are fitted to all the available. 
The parameters contain information on the fragmentation functions for the different quarks and 
the gluon and also a parametrisation of the non-perturbative contributions to the evolution. 
The value of the strong coupling constant obtained from this fit is 

asiMz) = 0.126 ± 0.007(exp) ± 0m6{theory) = 0.126 ± 0.009 . (7) 

The experimental error is the result of the combination in quadrature of the errors from the fit 
(0.0053), the uncertainties in the flavour composition of the enriched scaled energy distributions 
and the assumptions on the normalisation errors for those low-energy experiments where this 
error is not specified. The theoretical error is estimated by varying the factorisation and 
renormalisation scales. 

A possible extension of this analysis has been investigated by including the predicted dis- 
tribution measured at a centre-of-mass energy of 180 GeV (figure |]). Figure |10| shows the 
result of the fit to the scaled energy distributions at three centre-of-mass energies (29 GeV, 
91.2 GeV and 180 GeV). The fact that the variations with energy of the fragmentation func- 
tions is logarithmic makes the difference between the distributions at 180 GeV and 91.2 GeV 
smaller than that between 91.2 GeV and 29 GeV. This is accentuated by the fact that the 
flavour composition changes between 91.2 GeV and 180 GeV, in particular the percentage of b 
quarks diminishes when going to energies above the Z pole. Since the fragmentation function 
for b quarks is softer, this hardens the inclusive distribution at LEP II energies. 

The error in as{Mz) coming from the fit is not improved by including the distribution 
measured at 180 GeV. It was found, however, that with four times the predicted available 
statistics, a 10% improvement in this error could be obtained. The conclusion is that the 
analysis could serve as another consistency check of the predicted QCD scaling violations. 
Improvement in the error on as{Mz) may come from several sources. A better understanding of 
the flavour tagging algorithms used to measure the flavour-enriched distributions could improve 
the experimental systematic error. Progress on the theoretical side, for example the extension 
of the formalism to describe better the low-x region (see Section |5.3| ) could also be helpful. 



Another consistency check can be performed by using the measured flavour-enriched distri- 
butions at the Z peak and the scaling violation formalism to predict the fragmentation function 
in W decays. The fragmentation functions obtained from the fit to all data for the different 
quark flavours can be evolved to the mass of the W. Then the W scaled energy distribution 



can be predicted using the W decay branching ratios for each flavour. Figure |Tl| shows, in the 
continuous line, the prediction that results from this procedure. The points are the W scaled 
energy distribution as predicted by the PYTHIA Monte Carlo. 
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Figure 10: Result of the scaling violation fit to the distributions at centre-of-mass energies at 
29 GeV, 91.2 GeV and 180 GeV. 
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Figure 11: W fragmentation function predicted by the PYTHIA Monte Carlo (points), compared 
with the QCD prediction resulting from the analysis of scaling violations. 
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5.3 Small-x fragment at ionQ 



In the region of small values of the momentum fraction x, the behaviour of the fragmentation 
functions may be significantly affected by phenomena related to the coherence of soft gluon 
radiation (for a review of this subject see, for instance, Ref. [0). These effects are expected to 
result in a suppression of hadron production in the small-x (or soft) region, and to modify both 
the x-shape and the Q^-dependence of the inclusive single-particle spectrum. In particular, as a 
consequence of coherence, when the momentum fraction becomes small the gluon fragmentation 
function is expected to peak at a value dependent on the hard scale of the process, and be 
damped in the soft region. 

From the standpoint of perturbation theory, coherence effects show up as logarithmic cor- 
rections a|log™(l/x) {m < 2k — 2) to the splitting and coefficient functions which control 
the perturbative evaluation of the fragmentation functions. For example, the gluon splitting 
function Pgg{as,x) has the small-x behaviour {as = olsNc/t^) 

— — 2 — 3 

P{as,x)^ — -'^\og^x+'^\og^x + ... (8) 

XX 6X 

Small-x logarithms are present to all orders in as, and a systematic way to take coherence effects 
into account is to resum these logarithms to the leading accuracy, next-to-leading accuracy, and 
so on. 



The leading-log results were determined in Refs. [22, 47[, and can be best given in the 
moment space defined via the Mellin-Fourier transform 

jgg{as,uj)=[ dx x"^ Pgg{as,x) , (9) 

•/ 

and the analogous transform for any other function of x. In the moment space logarithmic 
terms appear as multiple poles at 0, and the summation of the leading contributions 



uj ^ ) is encompassed by the formula [^ 



7<?<?(«5, ^) = \ (^\/^^ + 8^5 - UJ^ . (10) 

The perturbative behaviour of this formula can be obtained by expanding it in the coupling 
0:5. The first terms of the expansion read as follows 

as a'i a% , , 

W«5,^)-^-2^ + 8-^ + ... , (11) 

where in the 0{as) and 0(a|) terms one may recognize the dominant part at small x of the 
standard one-loop and two-loop evolution kernels for the fragmentation functions (see [43] and 
references therein), whilst higher-order terms represent corrections due to coherent emission of 

^Author: F. Hautmann 
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soft gluons. An important feature which can be observed in Eq. (|TT|) is the alternating sign 
of the expansion. As a matter of fact, this feature extends to the whole series, and the net 
effect of resumming all the leading logarithms turns out to be a damping of the fragmentation 
function in the soft region with respect to the lowest-order prediction. 

The asymptotic properties of the resummed expression (0) are conversely determined by 
its behaviour near uj = 0. This is given by 



793 ("5, 
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Note that the all-order summation of the perturbative poles a^/uj'^'' ^ gives rise to a finite result 
at a; = 0, and introduces on the other hand the non-analytic behaviour in as of the square-root 
type. 



The summation of the next-to-leading contributions O (^ag/uj'^'^ has also been performed 
iD| . The explicit expression of the next-to-leading correction to Eq. ([1^) reads 



Yggi^^S.Uj) = 7gg + aS 
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where denotes the leading term (p!0D , and Nf is the number of flavours. Next-to-leading 
contributions do not alter the qualitative behaviour determined by the leading-order analysis, 
but provide a 0{y/as) correction to the position of the peak in the gluon fragmentation function. 



Phenomenological studies of the soft region of the single-particle spectrum have been carried 
out in Ref. |S0], on the basis of modified evolution equations which hold in the small-x regime. 
The central region of the spectrum, on the other hand, is known to be well described by second- 
order perturbation theory. It is therefore important to develop a procedure in which resummed 
contributions are consistently matched on to second-order perturbation theory, in order to get 
a uniform description of fragmentation over the whole phase space. 



6 Charged Particle Multiplicities^ 

The study of hadron multiplicity distributions in high energy collisions is an important topic 
in multiparticle dynamics and is generally undertaken as soon as a new energy domain be- 
comes accessible. It has been always considered a valuable tool to test our understanding 
of phenomenological approaches to multiparticle production and, in the framework of per- 
turbative QCD (MLLA) with assumption of Local Parton Hadron Duality (LPHD) the 

^Contributors: F. Fabbri and B. Poli (exp.), Yu.L. Dokshitzer and V.A. Khoze (th.) 
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average charged multiplicity < rich > and the second binomial moment of the distribution, 
R2 = ^"^""2''^ , are predicted to evolve with energy [20]. The measurement of the average 
charged multiplicity in heavy flavoured events is also of interest to perturbative QCD and a 
theoretical discussion on this particular topic is presented in this Section. 



6.1 Accompanying Multiplicity in Light and Heavy Quark Initiated 
Events 

Perturbative QCD approach predicts a suppression of soft gluon radiation off an energetic 
massive quark Q inside the forward cone of aperture Bq = Mq/Eq (Dead Cone)[^. This 
phenomenon is responsible for the "leading heavy particle effect" and, at the same time, induces 
essential differences in the structure of the accompanying radiation in light and heavy quark 
initiated jets. According to the LPHD concept [46], this should lead to corresponding differences 
in "companion" multiplicity and energy spectra of light hadrons. 

In particular, a solid QCD prediction is that the difference of companion mean multiplic- 
ities of hadrons, ANqi, from equal energy (hardness) heavy and light quark jets should be 
VT- independent [^, ^ (W is the energy available for soft particle production), up to power 



correction terms oc Mq/Wq. This constant is different for c and b quarks and depends on 
the type of light hadron under study (e.g., all charged, 7r°, etc). This is in a marked contrast 
with the prediction of the so called Naive Model based on the idea of reduction of the energy 
scale|^, Nqq(W) = Nqg{{l — {xq))W) , so that the difference of q- and Q-induced multiplicities 



grows with W proportional to N(W). 



The data[£6| for charged multiplicities in b- and c-quark events are in agreement with the 
energy independence of ANq£. As far as the the value of multiplicity differences is concerned, 
an expression for ANqe has been derived within the MLLA accuracy[^ assuming Mq ^ A: 



ANQe = Nqq{W) - Ngg{W) = -NggiV^MQ) [l + O (as(M2)) ] . (14) 

One usually consider the directly measurable quantity 

6Qi = ANq^ + Nd (15) 

where Nd is the average multiplicity due to the heavy quark decay. Quantitative QCD expec- 
tation for the difference of measured charged multiplicities ^qg that includes decay products of 
heavy hadrons, based on (|T^) was obtained in [^. For b quarks, which are only relevant for 



LEP2, the MLLA estimate 6m = 5.5 ± 0.8 exceeds the experimental value 2.90 ± 0.30. 



Recently an attempt has been made[^ to improve eq.(|Tj), the result of which modification 



agreed with the data ^^significantly better than the original MLLA prediction''^ (W.Metzger, 



56|). However, the very picture of accompanying multiplicity as induced by a single cascading 
gluon, implemented in [|^, is not applicable at the level of subleading O (as) effects (see, e.g. 
[46]). Therefore a reliable theoretical improvement of the QCD prediction for the absolute value 
of Nqi remains to be achieved. 
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DELPHI has recently measured ||59|| the number of vr^ in bb events to be close to that in all 
events. The same difference should be there at LEP2. 



6.2 Experimental 

The main limitations at LEP2 in this kind of study will come from the limited statistics and 
the relatively high contamination of events from other physical processes, absent or totally 
negligible at LEP I. Hadronic decays of W^W" pairs and highly radiative — > qq events 



are expected to constitute the dominant background. It was shown in subsection p.2| that 
this background can be reduced to a tolerable level, but the selection cuts needed, due to the 
particular nature of background events, will inevitably introduce a bias at both low and high 
multiplicities. 

The present study is based on the analysis of events generated with Pythia version 5.715 
at three different energies {^/s = 161, 175, 192 GeV), with hadronization parameters tuned to 
LEPl data [l^. The events were fully processed through the Opal detector simulation and 



reconstruction program chain, but the conclusions drawn here are believed to be practically 
the same for the other experiments. The statistics used was large compared to the most 
optimistic assumption on the integrated luminosity achievable at LEP2. Following the usual 
convention ||60|] , the charged multiplicity is defined as the total number of all promptly produced 



stable charged particles and those produced in the decays of particles with lifetimes shorter than 
3 ■ 10~^° sec. Non-radiative Z°/7 — »• qq events were selected following the criteria suggested in 
subsection B?^, in particular we used a combined "stage I" and Bj^ < 0.06 cut. A further 



background reduction was obtained by rejecting events with a Thrust value T < 0.8. The 
selection efficiency achieved for non-radiative events, defined as those with an Eisr < 1 GeV 
(see subsection p.2|) , was higher than 82% for all the considered energies. Due to detector 
acceptance and quality cuts, about 9% of the charged particles, on average, were lost in events 
surviving cuts while the predicted unbiased average charged multiplicity (< Uch >= 27.3 at 
^/s = 175 GeV) is about 14% higher than the observed one. Both those fractions were found 
to be practically energy independent. Approximately 33% of the events surviving cuts are 
radiative, namely with an Eisr > 1 GeV, but most of them have Eisr < 20 GeV. Background 
from W"'"W~ events never exceeds the 3% level. Residual background from other sources, like 
ZZ pairs, single W and single Z production, tau pairs and two-photon events was found to be 
negligible. 

The observed multiplicity distribution must be corrected for detector effects (acceptance 
and efficiency in track reconstruction, spurious tracks from photon conversions and particle 
interactions in the material, selection cuts) and for effects induced by the residual background. 



In figure |T2| we show the bias produced on the charged multiplicity distribution by the presence 
of residual W"'"W~ and radiative events as well as the bias produced by selection cuts. In fig- 
ure 



T^a we compare two normalized multiplicity distributions as they would appear in an ideal 
detector, namely without particle loss and interactions in the material, after event selection. 
In terms of real data, they would correspond to detector level corrected distributions including 
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residual initial state radiation (i.s.r.). The dotted distribution is relative to the pure qq sample 
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Figure 12: Bias from: a) residual W~^W ; b) radiative events; c) selection cuts. 

which survived cuts, the other distribution (histogram) contains also the residual contamina- 
tion from W"'"W~ events. As it can be seen from the bin-by-bin ratio shown in the bottom 
part, the bias due to this kind of background is relatively small. The estimated effect is 1.5% 



on < rich > and 0.6% on R2. In figure |T^-b a similar comparison is done to estimate the bias 
due to the presence of radiative events. Although this kind of contamination is relevant, the 
difference between the distribution containing residual radiative qq events (histogram) and the 
corresponding distribution for non-radiative events alone (dotted), is marginal. The effect is 
negligible on R2, while on < rich > is similar in size and in the opposite direction with respect to 
the one produced by W"'"W~ events. The bias introduced by selection cuts is important at high 
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multiplicity. This can be seen in figure ^^c where the normalized distribution for non-radiative 
qq events surviving the selection criteria, (histogram), and the normalized distribution for an 
unbiased qq sample, (dots), are compared. In this case the effect on < rich > and R2 was 
estimated to be 3.7% and 1.3%, respectively. 

Detector dependent corrections are usually carried out with unfolding matrix procedures 
and bin-by-bin coefficients [^]. In general matrices and coefficients are computed with the help 



of a very detailed detector Monte Carlo simulation in terms of material distribution, physi- 
cal processes which particles undergo when interacting in the material and detector response 
to particles traversing the active media. After subtraction of the estimated residual W"'"W~ 
contamination, using for example a bin-by-bin correction, and provided one has a reliable sim- 
ulation of the initial state radiation process, a global correction for particle loss due to detector 
effects is conceivable using a single unfolding matrix, computed from fully simulated events 
including i.s.r. The bias produced by event selection and residual i.s.r. can be corrected using 
bin-by-bin coefficients. 

Considering our estimated selection efficiencies and assuming cross sections and multiplicity 



distributions as predicted by Pythia, we show in figure |T3|-a,b the expected relative statistical 
uncertainties on < rich > and on R2 as a function of the integrated luminosity, at three different 
energies. The integrated luminosities expected at LEP2 are such that statistical uncertainties on 




Figure 13: Expected relative statistical uncertainties on < Uch > and i?2- 

these parameters should comfortably stay below the 1% level. A sensible estimate of the magni- 
tude of systematic uncertainties is difficult at this time. It will be most probably dominated by 
model dependent corrections needed to handle residual background and event selection biases. 
It is hard to believe it will be smaller than at LEPl (1 — 2% on < Uch >), but an uncertainty 
of a factor two higher may not be out of reach. We have fit the average charged multiplicity 
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measured above the Upsilon threshold[^ using the most popular parametrisations[p^] : 

< rich >= a - ■ exp(7/ 

< rich >=b- s"" 

< rich >= a + b ■ Ins + c ■ In^ s 

where a, b and c are free parameters. The predicted values extrapolated to LEP2 energies, 
however, only differ by a 3% or less and it will be not obvious to disentangle among these 
models. 

We also investigated the possibility to measure the average charged multiplicity in heavy- 
quark initiated events. Vertex-tagging methods have been shown to be very effective to select 
b-quark samples of high purity, and it is well known that secondary vertices with a relatively 
high associated multiplicity are likely to be produced in this kind of events. In the present study 
a method recently applied at LEP1|Q was used to analyse events simulated at -^i = 175 GeV. 
The method relies on the fact that independent samples of events with a different flavour 




Figure 14: a) Flavour composition vs. decay length significance; b) < rich > for the unbiased 
event hemisphere vs. b-quark purity. 

composition can be selected by requiring events to have (at least) one secondary vertex with 
a certain decay length significance, defined as the decay length divided by its error. The 
fraction of a given flavour in the sample can be evaluated, as a function of this variable, 
from fully simulated and reconstructed events, figure 0-a. In general to a high value of the 
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decay length significance corresponds a high probabihty to tag a b-quark event while at low 
significances the samples are predominantly populated by light-flavour events. To minimise the 
bias on multiplicity introduced by vertex tagging requirements, each event is divided in two 
hemispheres by a plane perpendicular to the thrust axis and the multiplicity is measured only 
in the hemisphere opposite to the one containing a secondary vertex, (unbiased hemisphere). 
The average charged multiplicity for pure samples of b-, c- and light-quark events is computed 
from a simultaneous fit to the corrected average multiplicity of samples selected with different 
decay length significance, i.e. with different flavour composition, figure |T^b. More details 



about the experimental procedure can be found in[34 



Due to extra selection requirements needed to insure the presence of secondary vertices, the 
original sample is reduced by a 30%. A b-tagging efficiency higher than 20% can be achieved 
in the highest purity bin. These values are similar to those obtained at LEPl. We studied 
the effects induced by the residual background and by the event selection cuts on the average 
charged multiplicity, measured in the unbiased hemisphere, as a function of the b-purity. Again 
we find that W~'"W~ and radiative events produce only a marginal effect while the bias produced 
by selection cuts is important. The unfolding matrices to correct for detector acceptance, 
efficiency and spurious tracks must be calculated for each bin of decay length significance. A 
bin-by-bin correction method could be used to unfold residual i.s.r. and selection cuts effects. 

In order to estimate the statistical precision attainable in this kind of measurement, we 
used the selection efficiencies found in this study and assumed the cross sections as well as 
the average multiplicity for different quark flavours, < Ug >, predicted by Pythia. The fltting 
procedure mentioned above was applied to a high statistics sample of unbiased qq events to 
estimate the uncertainty on < >. 



In flgure |T5l-a we show the expected relative statistical uncertainties on < rig > (q = b,c, light) 



as a function of the integrated luminosity. The difference in charged multiplicity between b- and 
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Figure 15: Expected relative statistical uncertainties on < Uq > and 6qi. 
light-quark events, 6bi, and between c- and light-quark events, 6ci, is shown, taking into account 
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correlations, in figure |l5|-b. One can see that these measurements will be largely dominated 
by statistical uncertainties, at least using this method of analysis. A measurement of 61,1 could 
probably be attempted, while a determination of 6ci seems to be precluded. 



7 Hadron Momentum Spectra as a Test of LLA QCD§ 

The shape of the momentum spectrum of hadrons produced in e'^e" collisions is successfully 
predicted in leading- log QCD (LLA). The LLA family of calculations together with the as- 
sumption of local parton hadron duality (LPHD) predict that soft gluons should interfere 
destructively due to their coherent emission (or angular ordering), and this gives rise to a 
'hump-backed' shape for the momentum distribution. At leading order, the distribution of 
ln(l/a;) should have a Gaussian form, and this shape is modified to be a Gaussian with higher 
moments if next-to-leading order terms are calculated. The position of the peak of the distri- 
bution ln(l/xo) is predicted in terms of the centre-of-mass energy, and therefore the evolution 
of the peak position with energy is well defined. This has been measured at centre-of-mass 
energies between 14 and 91 GeV, and the results are in agreement with a theoretical prediction 
that includes the effects of coherent, soft gluons. Models for hadron production based on phase 
space alone, or incoherent parton branchings predict a peak variation with energy that is twice 
as rapid and which is not supported by the data [62] |^ . 



The increase of centre-of-mass energy afforded by the energy upgrade LEP2, allows the 
hadron ln(l/x) distribution to be measured in a new energy regime and provides the opportunity 
to further test the evolution of ln(l/xo) from low energies. The energy increase is of the order 
of a factor two, so this represents a substantial 'lever arm' when compared to the existing data. 
In order to be able to challenge the predictions, the peak ln(l/xo) should be measured to a 
precision of less than about 0.1 unit of ln(l/xo). 

The detailed shape of the ln(l/x) is predicted in terms of a small (typically three) number of 
parameters, and an energy evolution. These parameters have been fixed by fitting to the LEPl 
data 1^ , and they can be used to predict the form of the data at higher energies. Clearly 
such a prediction of the shape of the ln(l/x) distribution constitutes an important potential 
measure of the success of the LLA approach to QCD calculations. 

The LLA approach has been extended to predict the momentum distribution of pairs of glu- 
ons which has commonly been presented in terms of the two particle correlation . This distri- 
bution has been measured at LEPl |6^[@ where it was found that the data were quahtatively, 
but not quantitatively described by next-to-leading order predictions. It was subsequently 
shown that a satisfactory description of the data was possible if next-to-next-to-leading 
order terms with coefficients of order unity were added to the prediction. The study of the 
energy dependence of the two particle correlation is interesting as it is predicted solely in terms 
of a single free parameter and the energy scale. 



Author: L.A. del Pozo 
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7.1 Monte Carlo Studies at 175 GeV 



Monte Carlo events generated at a centre-of-mass energy of 175 GeV were used to study the 
likely precision and limitations of an analysis of hadron momentum spectra at LEP2. Events 
were generated using Pythia for the processes e"'"e~ — >• W"'"W~ qqqq and e+e~ Z°/7 qq, 
and were subsequently passed through the Opal detector simulation program. 

In contrast to analyses at LEPl energies, the chief experimental problems are event statis- 
tics, and backgrounds in the event sample due to W+W~ qqqq events and events with a 
large amount of energy radiated by the initial fermions. The efficient selection of a clean sample 
of 2^/7 qq events with propagator energies close to the centre-of-mass energy has been 



extensively studied in subsection The present analysis uses the stage I cuts described there 
together with a stage II cut of > 300GeV^ and Bn < 0.05. In total about 9000 e+e" Z'^f-f 
events and about 130 e"'"e~ — >■ W^W~ events would be selected assuming the nominal LEP2 
luminosity of 500 pb~^ and standard model cross sections. About 72 % of the selected events 
have initial state radiation amounting to less than 2 GeV, and less than 2 % of the events have 
radiation in excess of 60 GeV. 



7.2 In {1/x) Distributions at 175 GeV 

The expected distribution of ln(l/x) is shown in figure [T6| (a) (and figure [16| (b) with a logarith- 
mic vertical scale) for all events that pass the selection cuts. The statistical errors on the points 
correspond to a luminosity of 500 pb~^. The contribution of the W"'"W~ background events is 
shown as the shaded region. The background is concentrated in the region around the peak 
of the ln(l/x) distribution, varies smoothly, and is small compared to the level of the signal 
events - the signal to background ratio is almost 100:1. In practice, this background could be 
corrected for by a multiplicative correction factor, the background could be subtracted directly 
or a more complex matrix correction procedure could be applied. 

Typical corrections for detector acceptance and resolution are shown in figure |16| (c). There 
is a correction of about 10 % to the overall level of the distribution which is basically fiat in 
the region around the peak of the ln(l/x) distribution. Uncertainties in the determination of 
detector corrections are therefore unlikely to have a large effect on the position of the peak, and 
there is no evidence that a serious bias has been introduced into the ln(l/x) distribution by the 
event selection. Figure |T^ (d) shows the ratio of the ln(l/a;) distributions for events passing the 
selection cuts that did not radiate and those that radiated a photon of more than 2 GeV. This 
illustrates the component of the detector correction that accounts for initial state radiation. 
The bias introduced by the initial state radiation is most severe for low values of ln(l/a;) but 
is fairly uniform around the area of the peak. It is not expected that the event selection and 
detector corrections will seriously bias the measurement of the position of the peak. 
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7.3 Determination of Peak Position 

The peak position may be determined by fitting a Gaussian to the W+W^-background sub- 
tracted ln(l/a;) distribution. The statistical error on the peak position is about 0.025 if data 
corresponding to 100 pb~^ are fitted, and this decreases to 0.011 when the expected 500 pb^^ 
data sample is analysed. A Gaussian function is only valid for the region close to the peak 
and is less successful at describing the shape of the distribution far away from it. This leads 
to a variation of the fitted peak position as data points far from the peak are included in the 
fit. Varying the fit range such that a reasonable is still obtained for the fit results in an 
uncertainty in the peak position of about 0.02. 

A systematic error due to uncertainties in the level of the W^W~ backgrounds has been 
estimated by varying the amount of the background subtracted by ±100 %. The fitted peak 
position changes by less then 0.01 in all cases. It is not expected that there is a large uncertainty 
due to the details of the shape of the W+W" background. The process W='= — > qq is very 
closely related to Z° — > qq which has been very well understood thanks the the LEPl data. 
In conclusion, it is expected that the position of the peak of the ln(l/a;) distribution may be 
measured with the data recorded at LEP2 to a sufficient precision in order to be able to test 
the LLA predictions. 

7.4 Detailed Shape of ln(l/a;) Distribution 

The expected statistical errors on the points of the ln(l/a;) distribution are small compared to 
the bin-to-bin variations of the distribution - there is no apparent scatter of the data points. 
This indicates that the data will most hkely be of sufficient precision to allow a detailed com- 
parison with the shape predicted by theoretical calculations with parameters fitted to LEPl 
data. The comments regarding the influence of acceptance, initial state radiation and W^W~ 
background corrections on the peak position also apply to the shape of the distribution - they 
are not expected to pose a major problem. If these systematic effects turn out to be trou- 
blesome, there is the still the prospect that a reasonable measurement of the ratio of ln(l/a;) 
distributions at LEPl and LEP2 energies might be made in which many systematic effects may 
cancel. 

It might be expected that the description of data by LLA predictions is more successful at 
higher energies as the LPHD assumption is more justifled. This is supported in part by Monte 
Carlo studies that indicate that the differences between hadrons and partons are much reduced 
at LEP2 energies. 
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7.5 Two Particle Correlation 



The two particle correlation at LEP2 energies has been studied in the same way as the single 
particle ln(l/a;) distribution. If the correlation distribution is computed along lines in the 
ln(l/a;i) — ln(l/x2) plane as in reference then the statistical error on each point would 
be of the order of 0.02 for the full 500 pb~^ data ample. This should be compared to 0.005 
achieved in reference with about 21 pb~^ of LEPl data. Preliminary studies indicate that 



corrections for acceptance, resolution and initial state radiation will be small as anticipated for 
this distribution. As for the LEPl analysis, it is also expected that other systematic effects 
such as the W"''W~ background might also cancel when the normalized correlation distribution 
is calculated. 

With the luminosity currently expected from LEP2 it is expected that any measurement 
of the two particle correlation would be statistics limited. There is however the hope that if 
the entire data sample is analysed, the possibility exists to test the energy evolution of the 
predicted correlation distribution in a meaningful way. In particular it can be tested whether 
the distribution at higher energies may be fitted by a prediction with the coefficients of the 
next-to-next-to- leading order terms fitted to LEPl data. Such a prediction with coefficients 
fitted to LEPl data is able to describe Pythia/Jetset events at both 91 and 175 GeV. Finally, 
the ratio of the two particle correlation at 91 and 175 GeV may be measured and compared to 
the theoretical prediction with the advantage that uncomputed higher order terms may cancel 
to some extent in the ratio. 



7.6 Summary 

In summary, measurements of hadron momentum spectra offer the possibility to make detailed 
tests of LLA QCD predictions, particularly in terms of their energy evolution. The peak 
position of the ln(l/x) distribution may be measured accurately with only a small amount 
of data allowing a powerful test of the extrapolation from lower energies. It should also be 
possible to determine the detailed shape of this distribution which will provide a stringent test 
of the energy evolution of predictions from LEPl energies. Meaningful measurements of the 
two particle correlation will probably have to wait for the full 500 pb~^ of luminosity to be 
dehvered by LEP2. 
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Figure 16: Distributions of ln(l/x) for 500 pb^^ of events passing selection cuts. Background 
from W^W~ events is shown as the shaded areas. Figures (c) and (d) show the the typical 
acceptance and initial state radiation corrections that might be expected. 
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